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We present a new symplectic structure and a hereditary recursion operator for the Sasa- 
Satsuma equation which is widely used in nonlinear optics. Using an integro-differential sub- 
stitution relating this equation to a third-order symmetry flow of the complex sine-Gordon II 
equation enabled us to find a hereditary recursion operator and higher Hamiltonian structures 
for the latter equation. 

We also show that both the Sasa-Satsuma equation and the third-order symmetry flow 
for the complex sine-Gordon II equation are bi-Hamiltonian systems, and we construct several 
hierarchies of local and nonlocal symmetries for these systems. 



I Introduction 

Finding a recursion operator for a system of PDEs is of paramount importance, as the whole in- 
tegrable hierarchy for the system in question is then readily generated by the repeated application 
of the recursion operator to a suitably chosen seed symmetry, see e.g. [7, 32, 11] and references 
therein. Moreover, using formal adjoint of the recursion operator enables us to produce infinitely 
many conserved quantities for our hierarchy (see e.g. [7, 11] for further details). 
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The recursion operators of multicomponent systems often have a richer structure of nonlocalities 
than in the one-component case, which makes such operators more difficult to find. For instance, there 
are two-component integrable generalizations, Eqs.(2) and (3), of the sine-Gordon equation and of the 
mKdV equation respectively. No recursion operator was found for Eqs.(2) and (3) so far even though 
the corresponding Lax pairs [36, 17, 38] and bilinear representations [17, 20, 18] are well known. The 
goal of the present paper is to fill this gap. Namely, we find and study the recursion operators for 
the so-called complex sine-Gordon II equation (2) and for the Sasa-Satsuma equation (3). 

The Sasa-Satsuma equation [36] (see also [46, 44]) for a complex function U has the form 

U t = U xxx + 6UUU x + 3U(UU) x . (1) 

Here and below the bar refers to the complex conjugate. 

It is natural to refer to this equation as to the complex mKdV II as (1) is one of the two inte- 
grable complexifications of the famous modified KdV equation, the other complexification (complex 
mKdV I) being simply U t = U xxx + 6UUU X . 

In turn, the complex sine-Gordon II equation [17, 38] is a hyperbolic PDE for a complex function 
ip of the form 

= i>Mv. + ( 2 ^ + C )(^ + C )hl>, 
ipip + c 

where c and k are arbitrary constants. The usual sine-Gordon equation <p xy = c 2 /csin(0) is recovered 
upon setting ip = ip = csin(0/4). 

The Sasa-Satsuma and the complex sine-Gordon II equation are of considerable interest for ap- 
plications. The Sasa-Satsuma equation is widely used in nonlinear optics, see e.g. [24, 34] and 
references therein, because the integrable cases of the so-called higher order nonlinear Schrodinger 
equation [25] describing the propagation of short pulses in optical fibers are related through a gauge 
transformation either to the Sasa-Satsuma equation [36] or to the so-called Hirota equation [19]. 

The complex sine-Gordon II equation, along with the Pohlmeyer-Lund-Regge model [33, 26, 
16, 30] also known as the complex sine-Gordon I, define integrable perturbations of conformal field 
theories [12, 2, 8, 3]; see e.g. [35, 37] for other applications. Moreover, the complex sine-Gordon I 
and II are the only equations for one complex field in the plane for which the (multi)vortex solutions 
are found in closed form [31, 5, 6]. 

Following [21] we set u = tp and v = ip and write the complex sine-Gordon II equation along with 
its complex conjugate as a system for u and v. 

u xy = — — - + (2-uf + c)(uv + c)ku, v xy = — + (2uv + c)(uv + c)kv. (2) 
uv — c uv + c 



Likewise, upon setting p = U and q = U in the Sasa-Satsuma equation, proceeding in the same 
fashion as above and writing out {pq) x as pq x + p x q, we obtain 

Pt = Pxxx + QpqPx + 3p 2 q x , qt = q X xx + 9pqq x + 3q 2 p x - (3) 

From now on we shall treat u, v, p and q as independent variables that can be real or complex 
and consider systems (2) and (3) that are more general than the original complex sine-Gordon II 
and Sasa-Satsuma (complex mKdV II) equations which can be recovered under the reductions v — u 
and q = p respectively. In what follows we shall refer to (2) as to the complex sine- Gordon II system 
and to (3) as to the Sasa-Satsuma system. 

Let us briefly address the relationship of the Sasa-Satsuma system (3) with other integrable 
systems. First of all, Eq.(3) can be obtained [46] as a reduction of the four-component Yajima- 
Oikawa system (8) from [46]. 

On the other hand, consider the vector modified KdV equation 

V t = V xxx + (V, V)V X + (V, V X )V (4) 

studied by Svinolupov and Sokolov [42] (see also [43, 45]). Here V = (V 1 , . . . , V n ) T is an n-component 
vector, (•, •) stands for the usual Euclidean scalar product of two vectors, and the superscript T here 
and below indicates the transposed matrix. A simple linear change of variables V 1 = y/3(p + g)/\/2, 
V 2 = iV3(p - q)/y/2 takes (3) into (4) with n = 2. 

Setting r = \^3p and s = V3q turns (3) into the system 



< XXX 

+ r(3sr x + rs x ), 
s t s xxx ~l" s (3rs :r + st-j,), 
studied by Foursov [15], who found a Hamiltonian structure for (5) of the form 



(5) 



P = 



( -\rD^or D x + \rD- l o S \ 



(6) 



\D x + \sD~ 1 or -IsD-^os J 

where D x is the operator of total x-derivative, see e.g. [7, 11, 32] for details and for the background 
on the recursion operators, Hamiltonian and symplectic structures. The corresponding Hamiltonian 
density is [15] (2/3)r 2 s 2 — r x s x . This immediately yields a Hamiltonian structure (8) for (3). Note 
that in [14, 15] there seems to be a misprint in P, and in (6) we corrected this misprint. 

In [14, p. 89] Foursov claims to have found some skew-symmetric operators that are likely to 
provide higher Hamiltonian structures for (5) but he failed to verify that these operators are indeed 



Hamiltonian. The explicit form of these operators was not presented in [15] or [14], so apparently it 
was never proved in the literature that (5) (and hence (3)) are bi-Hamiltonian systems. We establish 
the bi-Hamiltonian nature of (3) in Theorem 1 below. 

Now turn to the complex sine-Gordon II system (2). Recall that Eq.(2) can be obtained (see e.g. 
[17, 21, 28]) as the Euler-Lagrange equation for the functional S — f L dx dy, where 

1 U X Vy + UyV X 

L = - - h k(uv + cjuv. 

2 uv + c v ' 

A few conservation laws and generalized symmetries for (2) can be readily found e.g. using computer 
algebra [21, 28]. In particular [21], Eq.(2) is compatible with 

_ _ 3uv x u xx _ 9u 2 x v x 3u 2 v 2 u x 

U>t t^xxx , , + 



UV + C uv + c (uv + c) 2 ' 

(7) 

_ 3vu x v xx 9vlu x 3v 2 u 2 x v x 

Vt — v xxx + . . . 

uv + c uv + c [uv + cy 
The compatibility here means that the flow (7) commutes with the nonlocal flow (14) associated 

with (2) or, equivalently, the right-hand sides of (7) constitute the characteristic of a third-order 
generalized symmetry for (2), see e.g. [32] for general background on symmetries. 

The rest of paper is organized as follows. In section II we present symplectic structure and re- 
cursion operator for the Sasa-Satsuma system (3) and show that (3) is a bi-Hamiltonian system. In 
section III we employ a nonlocal change of variables relating systems (7) and (3) in order to construct 
recursion operator, Hamiltonian and symplectic structures for (2) and (7) from those of (3), and we 
show that (7) is a bi-Hamiltonian system. Finally, in section IV we discuss the hierarchies of local 
and nonlocal symmetries for (2), (3) and (7). 

II Recursion operator and symplectic structure 
for the Sasa-Satsuma system 

A straightforward but tedious computation proves the following assertion. 
Theorem 1 The Sasa-Satsuma system (3) possesses a Hamiltonian structure 

' -P D x l °P D x +pD~ 1 oq 
yD x + qD- 1 op - q D- l oq 

' 3pD~ l op D x + 5pD- l oq 
y D x + 5qD- 1 op SqD^oq 
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V 

a symplectic structure 

J = 



(8) 



(9) 



and a hereditary recursion operator 1Z = V o J that can be written as 

( D\ + Qpq + q x D~ x op + bD x x op + 3p x D- x ob 2p 2 - 2z 1 pD x 1 oq- 3p x D x 1 o a ^ 
^ 2q 2 - 2z 2 qD~ 1 op + 3q x D~ 1 ob D 2 X + Qpq + p x D~ x oq + aD' 1 oq- 3q x D~ 1 o a 



ft 



(10) 



where a = p x + 2z\q, b = q x + 2z 2 p; z\ = D x 1 {p 2 ) and z 2 = D x 1 (q 2 ) are potentials for the following 
conservation laws of (3): 

D t {p 2 ) = D x (2pp xx -p 2 x + 6p 3 q), D t {q 2 ) = D x (2qq xx - q 2 x + Qpq 2 ). 

Hence (3) has an infinite hierarchy of compatible Hamiltonian structures Vk = Tl k ° V , k = 
0,1,2,..., Vo = V , an infinite hierarchy of symplectic structures J k — J o 1Z k , k = 0,1, 2,..., 
and an infinite hierarchy of commuting symmetries of the form /Q = 1Z % (JCq), i = 0, 1, 2, . . where 

= (Px,qx) T - 

The Sasa-Satsuma system (3) is bi- Hamiltonian with respect to Vq and V\. 







( 5Hi/5p \ 




' 5Ho/Sp N 








= v 1 








y 5Hi/5q J 




^ sn /Sq j 



where Hi = J Hidx, i = 0, 1, H = pq, H x = 2p 2 q 2 - p x q x . 

Here 5/5p and S/Sq denote variational derivatives with respect to p and q. 

Note that the Hamiltonian structure V above is nothing but the second Hamiltonian structure 
of the AKNS system, see e.g. [27]. This Hamiltonian structure (more precisely, its counterpart (6) 
for Eq.(5)) has already appeared in [14, 15]. 

The symmetries /Q commute because 1Z is hereditary. Applying 1Z to an obvious symmetry 

= (Px, q x ) T of (3) yields the symmetry /Ci = K(Ko) = (p X xx + QpqPx + 3p 2 q x , q X xx + Qpqqx + 3q 2 p x ) T , 
i.e., the right-hand side of (3). In turn, TZ(ICi) is a local fifth-order symmetry for (3). We guess that 
/Q are local for all natural i but so far we were unable to provide a rigorous proof of this. 

Unlike the overwhelming majority of the hitherto known recursion operators, see e.g. discussion in 
[39] and references therein, the nonlocal variables appear explicitly in the coefficients of 1Z. Perhaps 
this is the very reason why 1Z was not found earlier. The nonlocal variables in the coefficients of 
1Z are abelian pseudopotentials as in [23] and unlike e.g. the nonlocalities in the recursion operator 
discovered by Karasu et al. [22] and later rewritten in [40]: the nonlocalities in the operator from 
[22, 40] are nonabelian pseudopotentials. 



Ill Recursion operator 

for the complex sine-Gordon II system 



There is a well-known (see e.g. [1] for discussion and references) transformation z = ^2/3g x relating 
the symmetry flow g t = g xxx + g x i /2 of the sine-Gordon equation g xy = sin(g) and the mKdV equation 
z t = z xxx + z 2z x- Moreover, there exists [1] a nonlocal generalization of this transformation that sends 
the third-order symmetry flow 



_ _ Q uv x u xx {-uvu x - cu x + u 2 v x )u x v, 

— Uxxx <J . + <J 



uv + c (uv + c) 2 



X 

1 



_ vu x v xx (-v 2 u x + uvv x + cv x )u x v x 

v t — v xxx "J , <J _\ 2 



uv + c (uv + c) 

of the complex sine-Gordon I equation (see e.g. [16, 21, 28]) 

u xv = h (uv + c)ku, v xv = — + (uv + cjkv, (11) 

uv + c uv + c 

into the two-component generalization of the mKdV equation 

Pt = Vxxx + Qpqpx, qt = qxxx + Qpqqx 

that belongs to the hierarchy of the well-known AKNS system 

Pt = Pxx + p 2 q, qt = -q xx - q 2 p- ( 12 ) 

Note that this nonlocal transformation also sends [1] the second-order symmetry flow 

_ uu x v x vu x v x 
u t — u xx — * ; — , Vt — —v xx + z- 



uv + c uv + c 

of (11) into the AKNS system (12). 

It turns out that upon a suitable redefinition of nonlocal variables the nonlocal transformation 

in question in combination with a suitable rescaling of dependent variables p and q also sends the 

third-order symmetry flow (7) of the complex sine-Gordon II equation (2) into the Sasa-Satsuma 

system (3). Namely, we have the following result. 

Theorem 2 The substitution 

iv / 2Mxexp(-iwi) iV2v x exp(^Wi) 
P= . — — , q = . = — , 

y/UV + C y/UV + C 

where i = \J — \, W\ = D~ 1 (p 1 ) is a potential for the conservation law D t (p 1 ) = D x (a\) of (7), 

uv x — vu x uv xxx — vu xxx (3uv + 2c)(v x u xx — u x v xx ) 
Pi = ; > °\ = ; 1- 



uv + c ' uv + c (uv + c) ; 

12uv + 11c) (vu x — uv x )u x v x 



x 

i 



(uv + c) 3 

takes Eq.(7) into the Sasa-Satsuma system (3). 



Here and below D x and D t stand for the total x- and t-derivatives as defined e.g. in [32]. 
Remark 1 Let 



P 



iv / 2Mxexp(-|wi) 



-, q 



iV2v x exp(^w 1 ) 



^/uv + c ' \Juv + c 

where w\ = D~ 1 (p\) is now a potential for the conservation law D y (pi) = D x (6\) of (2), p\ is given 
above and 

vu y — UVy 



0i 



uv + c 



Then for k = we have 



D„ 



0, D y (q) = 0, 



where D y stands for the total y-derivative. In other words, if k = then p and q provide nonlocal 
y -integrals for (2), and D x Yn{p) and D x \n(q) are local y -integrals of (2). Therefore, the complex 
sine-Gordon II system (2) for k = is Liouvillean and C-integrable, see [9, 10] for the construction 
of symmetries of such systems using the integrals thereof. Using the above local integrals enables us 
to find the general solution for (2) with k = along the lines of [47]. 

Passing from p and q to u and v yields from 1Z a recursion operator !tH for (7). The operator 
is hereditary because so is 1Z, see [13]. It is easily seen that 9^ is a recursion operator for (2) as well. 
Upon removing an inessential overall constant factor in we obtain the following result: 



Theorem 3 The complex sine- Gordon II system (2) has a hereditary recursion operator 



( p 2 2uv x t uv xx - 12u x v x t cu x v x -2uu xx + Qu 2 x i Au 2 u x v x \ 



+ 



V 



uv + c uv + c 

+ 



+ 



(uv + c)' 



uv + c 



+ 



(uv + c) 2 



2vv xx + 6v 2 | 4v 2 u x v x ^ 2 2vu x ^vu xx — \2a r r, cu ,.r 



uv + c 



(uv + c) 5 



uv + c 



uv + c 



+ 



(uv + c) 2 / 



+ 



( E QlD- 1 o 7a>1 £ Q\ D x X ° 7a>2 X 

a=l a=l 



E QlPx 1 ° 7a,i EWov 

a=l 



\ a=l 



where w 1 is as in Remark 1; y\ = D x 1 ((-u 2 exp(— Wi))/ (uv + c)) and y 2 = D x l ((v 2 exp(wi))/(uv + c)) 
are potentials for the nonlocal conservation laws D y (Q) = D x (£i), i = 1,2, of (2), 

. ulexp(-wi) ^ , 9 . . . . 
Ci = ... , . , £i = ku (uv 4~ c) exp( — tfi), 



Mf + c 

^exp(wi) 
Mf + c 



, £ 2 = kv 2 (uv + c) exp(wi); 



Qi 



yiu 

y w x exp(-wi) - vy x J 



\ ( 

, Q 2 = 



( u ^ 



\- v J 



, Q 4 = c 



/ 



V 



-U a 



^exp(wi) - My 2 
2uu x v x 



uv + c 



Aexp(w 1 )y 1 v x - 4wy 2 wi 



2vu x v x . . 

■ 4 exp(wi)y 2 Mx + 4vy 2 yi / 

Mt> + c / 



5 



Q 5 = 



V 



are symmetries for (2), and 

' exv(w 1 )vvl + cy 2 v x exp(w 1 )v xx 2exp(w 1 )vu x v x + cy 2 u a 



7i 



7 2 = 



(m1> + c) 2 Wf + C (wt> + c) 2 

exp(-M>i)M x;r 2exp(-iui)uu x 7; x + ^1^ exp(-M>i)MM 2 + q/iM x 



uv + c 



(uv + c) 2 



(m1> + c) 5 



7 3 



+ 



2(w x + 2exp(-wi)ttv?/2 + 2exp(-M>i)cM 2 )M xx + (2vw x + uv x )v x 



uv + c (uv + c) 2 

2v x u x (-v 2 u x + 2cv x ) 2v x (-u x v x + 4exp(-w 1 )uy 2 u x + 2 exp(w 1 )vy 1 v x + 2cy 2 yi) 



uv + c 



(uv + c) 3 (uv + c) 2 

_ 2(w^ + 2exp(w l )uvy 1 + 2 exp(wi)cy 1 )v xx - (2uv x + vu x )u a 

(uv + c) 2 



2v x u x (-u 2 v x + 2cu x ) + 2M x (-w x t; x + 4exp(wi)tn/it; x + 2 exy(-wi)uy 2 u x + 2q/ 2 Mi) \ 



(mi; + c) 3 



(uv + c) 5 



74 = 



(■uf + c) 2 ' (wf + c) 2 / 
are cosymmetries for (2). 



7 5 



M 7 



vu. 



(uv + c) (mi; + c) 2 ' (mi; + c) (mi> + c) 5 



Here Q l a and 7 a j denote i th component of Q a and j th component of 7 a respectively, i.e., 



) 7a — (7a, 1? 7a, 2 



Note that if we use the tensorial notation (see e.g. [28]), we can rewrite the operator 

/ 5 5 \ 



a=l a=l 



1 EWov EWov . 

\ a=l a=l / 



in a more concise form, namely Yl Qa ® -D x 1 ° 7 Q - 

a=l 
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Remark 2 Eq.(7) has a recursion operator of precisely the same form as the fR given above, but 
in this case the nonlocal variables should be defined in a slightly different way: W\ should be as in 
Theorem 2, and yi = D~ x ((w 2 exp(— wi))/(uv + c)), y 2 = D~ l ((v^.exp(wi))/(uv + c)) should now 
be potentials for the nonlocal conservation laws D t (Q) = D x (xi), i = 1,2, of (7), where Q are as in 
Theorem 1 and 



Xi = exp(-wi) 



X2 = expfwi, 



uv + c 

2v x V X xx 



U , 



2uu 2 x v xx 



V. 



uv + c 

2 



uv + c) 2 (uv + c) 2 
2vu x v x v xx , 



2(6uv + 7c]u\v x 3u 2 ulv 2 



2vv\u xx 



{uv + c) 3 (uv + c) 3 

2 "' 2 "' 2 2(6uv + 7c)u x v 



3v u x v x 



uv + c uv + c (uv + c) 2 (uv + c) 2 



[uv + c) 



(uv + c) 3 



Recall that Eq.(2) possesses [28, 29] a local symplectic structure 

/ 1 „ uv r 

3 = 







-D x 



uv + c 
vu x 



-D,, 



(uv + c) ; 







\ 



/ 



\ uv + c x ("Uf + c) 2 

i.e., Eq.(2) can be written in the form %u y = k(2uv + c)(v,u) T = 5H/d~u, where u = (u,v) T and 
H = J k(uv + c)uvdx. The symplectic structure 3 is readily seen to be shared by Eq.(7). 
Inverting 3 yields a nonlocal Hamiltonian structure for (7) of the form 

^ exp(— wi/2)y/uv + co D~ x o exp(wi/2)y/uv + c 

y exp(wi/2)^uv + co D~ l o exp(— wi/2)y/uv + c 

where w\ should be interpreted as in Theorem 2. It can be shown that if we go from the (u, v)- to the 



(p, g)-variables using the transformation from Theorem 2 then we obtain from *P the Hamiltonian 
structure (8) for the Sasa-Satsuma system up to an inessential overall constant factor. 

The higher Hamiltonian (resp. symplectic) structures for (7) are given by the formulas = 
9t n o <p (resp. Zn = 3 ° 9^ n ), where n = 1, 2, 3, . . .. The operator !tH is hereditary, and hence all of 
these structures are compatible. However, the structures tyi and 3i are already very cumbersome, 
so we do not display them here. Nevertheless, it is readily checked that the following assertion holds. 

Theorem 4 Eq.(7) is a bi- Hamiltonian system: 







/ 


8f)i/5u 


\ 




( 


5S) /5u 


\ 




















V 


5$ii/5v 


J 




\ 


5fi /5v 


J 



where 9)i = f Hidx, i — 0, 1, 



uv + c 



U^xx^xx VU x V x Ux X ~t~ UU X V X V XX ^^x^x 

UV + c (uv + c) 2 



uvu 2 x v 2 
(uv + c)' 
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It is straightforward to verify that the action of fR on the obvious symmetry u x of (7) yields 
the right-hand side of (7), and 9\ 2 (u x ) is a fifth-order local generalized symmetry for (2). We guess 
that the repeated application of £H to u x yields a hierarchy of local generalized symmetries for (2) 
but we were not able to prove this in full generality so far, as presence of the nonlocalities yi in the 
coefficients of £H appears to render useless the hitherto known ways of proving locality for hierarchies 
of symmetries generated by the recursion operator, cf. e.g. [41] and references therein. Nevertheless, 
as is a recursion operator for (7), Proposition 2 from [39] tells us that the only nonlocalities that 
could possibly appear in the hierarchy of the symmetries Vi k (u x ), k = 1,2,..., are potentials of 
(possibly nonlocal) conservation laws for (2). 

IV Nonlocal symmetries 

First of all, we can consider Eq.(2) as a nonlocal symmetry flow of Eq.(7). We have already noticed 
above that Eq.(2) can be written as 



where Sj = f k(uv + c)uvdx. Acting by ^3 = Z 1 on both sides of (13) we can formally rewrite the 
complex sine-Gordon II system (2) in the evolutionary form u y = ty(5$)/5u) : that is, 



(13) 



u y = k exp(— w\/2)\/uv + c tu 2 , v y = keyL^{wi/2)\fmT^rc ui, 



(14) 



where u\ = D x 1 (y(2uv + c) exp(— w 1 /2)^uv + c) and uo 2 = D x 1 (u(2uv + c) exp(wi/2)y/uv + c) are 
potentials for the following nonlocal conservation laws of (7): 




2v 3 u 2 x v(20u 2 v 2 + 26cm; + 7c 2 )u x v : 




uv + c (uv + c) 2 




, w(20w V + 26cm; + 7c 2 )u x v x 2u 3 v 2 \ \ 

~ 2vu x 7 ^ " ~ 

[uv + c) z uv + c / / 

Eq.(14) is compatible with Eq.(7), i.e., the right-hand side of (14) divided by k, 




Q i = 



y exp(wi/2)^uv + c lo\ J 
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is a nonlocal symmetry for (7). Moreover, we have 9 I l(Q_ 1 ) = c 2 u X) i.e., the action of 9\ on Q_ 1 
gives the 'zeroth' (obvious) symmetry Uj. up to a constant factor. Thus, Eq.(14) (and hence the 
complex sine-Gordon II system (2)) can be considered as a first negative flow in the hierarchy of (7). 
Note that if we consider £H as a recursion operator for Eq.(2), we have 9K(u y ) = kc 2 u x , in perfect 
agreement with the above result. 

Eqs.(3) and (2) possess nonlocal symmetries of the form 

/ \ 



Qi 



^ 1x + Z 2 p ^ 



and 



Qi = 



yiu 



-z 2 q 



J 



u x exp{-w 1/ 



vyi 



J 



Zip 

( 



Q 2 = 



v x ex.p(w ± ) - uy 2 



\ 



\ 



J 



respectively, where Z\ and z 2 were defined in Theorem 1. Therefore we have four hierarchies of non- 
local symmetries: TZ l (Gj) for (3) and 9t l (Q 3 -) for (2), where in both cases j — 1,2 and % — 0, 1, 2, 

Two more hierarchies of nonlocal symmetries have the form TZ l (Go) fc> r (3) and 9i l (Q ) for (2), 



where i 



1,2, 



Go 



( " ) 

\- q J 



( 



u 



\ 



Two somewhat more 'usual' nonlocal hierarchies of master symmetries (the latter are, roughly 
speaking, time-dependent symmetries such that repeatedly commuting them with a suitable time- 
independent symmetry yields an infinite hierarchy of time-independent commuting symmetries for 
the system in question, see e.g. [7] and references therein for further details), TZ l (So) for (3) and 
9V 1 (So) for (7), where i — 1, 2, . . ., originate from the scaling symmetries for (3) and (7), 

( 3t(p xxx + 9pqp x + 3p 2 q x ) + xp x + p 

,2„ 



S = 



and 



y 3t(q xxx + 9pqq x + 3q p x ) + xq x + q 

3 UV x 1^ x x 



So = 



( 



\ 



3t u 



9u x v x 3u v x u x 
uv + c uv + c (uv + c) 2 



+ xu x 



\ 



3t v 



^DIlx'Vxx 



9v 2 x u x 3v 2 u 2 x v x 



uv + c uv + c (uv + c) 2 



XVn 



respectively. 

Finally, (7) has nonlocal symmetries 

^ \/uv + cexp(— Wi/2) 



G 1 = 



Go 



\ 



\fwv + c exp (wi 1 2) 
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obtained by differentiating Q_ x with respect to uj\ and lo 2 , but these symmetries are annihilated by 
and hence do not lead to new hierarchies of nonlocal symmetries. 

It would be interesting to find out whether the systems in question possess nonlocal symmetries 
that do not belong to the above hierarchies, which is the form of solutions invariant under the 
nonlocal symmetries and whether these solutions could have any applications in nonlinear optics. 

As a final remark, we note that because of the obvious symmetry of the complex sine-Gordon II 
system (2) under the interchange of x and y, all of the above results concerning the recursion opera- 
tor, Hamiltonian and symplectic structures, and hierarchies of symmetries for (2) remain valid if we 
replace all x-derivatives with y-derivatives and vice versa and swap the operators D x and D y . Interest- 
ingly enough, the recursion operator 9^ obtained from 9^ upon such an interchange proves to be inverse 
to 91 on symmetries of (2) up to a constant factor. More precisely, for any symmetry K of (2) we have 

&(K) = k 2 c 4 9T\K). 

Taking into account our earlier results we see that the 'basic' hierarchy of (2) can be represented by 
a diagram of the form 

and we guess that all symmetries presented at this diagram are local, i.e., they do not involve nonlocal 
variables. This is a fairly common situation for hyperbolic PDEs, see e.g. the discussion in [4]. 
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